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Abstract:  This  paper  describes  an  algorithm  using  a  discrimination-based  sensor 
Effectiveness  metric  for  sensor  assignment  in  multisensor  multitarget  tracking 
applications.  The  algorithm  uses  Interacting  Multiple  Model  Kalman  Filters  to  track 
airborne  targets  with  measurements  obtained  from  two  or  more  agile-beam  radar  systems. 
Each  radar  has  capacity  constraints  on  the  number  of  targets  it  can  observe  on  each  scan. 
For  each  scan  the  expected  discrimination  gain  is  computed  for  the  sensor  target  pairings. 
The  constrained  globally  optimum  assignment  of  sensors  to  targets  is  then  computed  and 
applied.  This  is  compared  to  a  fixed  assignment  schedule  in  simulation  testing.  We  find 
that  discrimination  based  assignment  improves  track  accuracy  as  measured  by  both  the 
root-mean- square  position  error  and  a  measure  of  the  total  covariance. 

1.  Introduction:  The  goal  of  sensor  management  systems  can  be  regarded  as  selecting 
sensors  and  sensor  dwells  to  increase  the  information  that  a  data  fusion  system  contains 
about  a  surveillance  region.  Discrimination  [Blahut]  is  related  to  the  notions  of 
information  and  entropy  in  probability  distributions.  A  sharply  peaked  distribution  has 
high  discrimination  gain  and  low  entropy.  Discrimination  measures  the  relative  increase 
in  information  between  two  probability  distributions.  In  light  of  this,  discrimination  gain 
is  a  natural  measure  for  use  in  sensor  management  systems.  Expected  discrimination  gain 
is  a  measure  of  sensor  effectiveness  that  has  been  used  in  a  wide  variety  of  model 
applications  including  multisensor  /  multitarget  assignment  problems  [Sch93], 
minimizing  error  correlation  between  close  targets  [Sch94],  and  single  and  multi  sensor 
detection  /  classification  problems  [Kas96b,  Kas97a]  and  joint  detection,  tracking  and 
identification  [Kas96c,  Kas97b].  It  compares  favorably  to  other  methods  [Mus96].  In  the 
approach  presented  here  this  entails  predicting  the  expected  discrimination  gain  for  each 
sensor  dwell,  determining  the  optimal  global  assignment  of  sensors  to  targets  given 
sensor  constraints  and  then  applying  this  solution. 

To  understand  how  expected  discrimination  gain  can  be  computed  in  Kalman  filter  based 
tracking  systems,  recall  that  Kalman  filters  maintain  estimates  of  both  the  target  state  and 
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the  state  covariance.  The  target  state  consists  of  the  target  position  and  velocity.  It  may 
also  include  higher  order  kinematics  terms  such  as  turn  rate  and  acceleration.  The 
covariance  is  a  measure  of  the  uncertainty  in  these  state  estimates.  In  addition  to  the  state 
and  covariance  estimates  computed  after  each  measurement  update,  the  Kalman  filter  also 
includes  a  mechanism  to  predict  the  state  at  future  times  based  on  assumptions  about  the 
target  motion  (the  "motion  model").  The  uncertainty  in  this  prediction  can  also  be 
computed.  Given  the  predicted  covariance,  the  expected  discrimination  gain  can  then  be 
evaluated.  A  surprising  feature  of  conventional  Kalman  filters  is  that  the  covariance 
estimate  depends  only  on  the  uncertainty  of  the  measurements,  the  state  uncertainty  and 
the  target  motion  model,  but  it  does  not  depend  directly  on  the  measurements  or  the  target 
state  itself.  (In  an  extended  Kalman  filter,  the  uncertainty  can  depend  on  the 
measurements  indirectly  through  the  linearization  process.)  Thus,  in  conventional 
Kalman  filters  the  uncertainty  and  the  expected  discrimination  gain  are  insensitive  to 
whether  a  target  is  maneuvering  or  not.  This  undesirable  feature  reduces  the  effectiveness 
of  sensor  management  systems  based  on  them. 

The  insensitivity  of  the  Kalman  filter  to  target  maneuvers  has  long  been  recognized  and  a 
wide  variety  of  methods  have  been  developed  to  recognize  maneuvers  and  change  the 
filter  accordingly.  Recently,  a  new  approach  called  the  Interacting  Multiple  Model 
Kalman  Filter  (IMMKF)  has  been  developed  that  seems  to  offer  superior  performance  to 
other  approaches  [Blom88,  Blom92].  The  IMMKF  uses  several  target  motion  models 
and  has  been  successfully  applied  to  large  multisensor  multitarget  applications  [Yed]. 

For  example,  it  may  use  one  model  for  straight  and  level  flight  and  different  models  for 
turns.  The  turning  models  can  use  a  fixed  turn  rate  [Li]  or  can  estimate  the  turn  rate  as  an 
additional  state  in  the  filter  [Kas96].  Older  approaches  to  filtering  switch  between 
models.  In  contrast  to  this,  the  IMMKF  always  maintains  all  of  the  models  and  blends 
their  outputs  with  weights  that  are  computed  probabilistically.  In  addition  to  the  state 
estimates  for  each  motion  model  the  IMMKF  maintains  an  estimate  of  the  probability  that 
the  target  is  moving  in  accordance  with  each  model.  With  the  IMMKF  the  prediction 
covariance  increases  during  target  maneuvers.  As  a  result  the  expected  discrimination 
gain  goes  up  when  a  sensor  is  directed  at  a  target  that  is  perceived  to  be  maneuvering. 
Another  factor  that  contributes  to  the  expected  discrimination  gain  with  the  IMMKF  is  its 
estimate  of  which  target  motion  model  is  active.  As  a  result,  additional  sensor  resource  is 
directed  toward  the  target  as  it  is  initiating  a  maneuver. 

The  natural  coupling  between  the  IMMKF  and  sensor  management  for  single  target 
applications  was  recognized  in  [Watson]  and  exploited  using  measures  of  effectiveness 
based  on  maintaining  a  desired  level  of  position  error  covariance.  Coordination  across 
multiple  sensors  is  desirable  in  multitarget/multisensor  applications.  Typically,  each 
sensor  can  only  observe  a  few  targets  on  each  sensing  cycle.  By  coordinating  sensing 
across  platforms,  sensing  effort  is  not  wasted  by  having  several  sensors  expend  all  of  their 
capacity  on  one  target  while  other  targets  go  unobserved. 

The  optimization  problem  considered  here  is  one  of  assigning  sensors  to  targets  in  such  a 
way  that  the  one  step  ahead  discrimination  gain  of  the  track  set  is  maximized.  To  do  this, 
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the  one  step  ahead  expected  discrimination  gain  Gsl  of  a  track  t  after  it  is  updated  with 
covariance  data  from  sensor  s  is  determined  for  all  pairs  t,s.  This  can  be  formulated  as  a 
Linear  Programming  problem  with  suitable  constraints.  Situations  where  more  than  one 
sensor  may  be  assigned  to  the  same  target  can  also  be  treated  [Sch93].  As  has  been 
previously  pointed  out  [Nash]  one  method  is  to  construct  pseudo  sensors  comprised  of 
combinations  of  the  basic  sensors.  This  allows  any  combination  of  sensors  to  be  assigned 
to  a  single  target  simply  by  considering  that  a  single  pseudo  sensor  has  been  assigned  to 
it.  The  number  of  "sensors"  is  thus  increased  from  S  to  2s  -  1.  The  problem  is  now  to 
make  an  assignment  of  these  pseudo  sensors  to  the  targets  in  an  optimal  way. 

There  are  however  constraints  on  this  assignment.  One  of  the  most  important  is  the 
maximum  tracking  capacity  of  a  sensor.  That  is,  given  a  specified  time  interval,  some 
sensors  can  scan  a  certain  volume  of  space  and  also  track  a  specified  number  of  targets. 
They  cannot  exceed  this  maximum  tracking  capacity  in  the  specified  time  period.  To 
examine  the  impact  of  coordination  between  the  sensors,  we  studied  three  cases.  In  the 
first  case,  each  sensor  simply  selects  targets  at  random  for  each  scan  cycle.  This  is  the 
worst  case  approach  in  that  it  assumes  that  the  sensing  is  completely  uncoordinated.  This 
can  be  improved  upon  by  assuming  that  the  sensing  is  coordinated  using  a  round-robin 
sensing  schedule.  The  result  of  these  approaches  is  compared  to  the  global  discrimination 
optimum.  To  establish  a  baseline  for  comparison,  all  of  the  scheduling  schemes  maintain 
central  tracks  updated  via  report  fusion  as  soon  as  measurements  are  available.  Although 
this  may  be  unrealistic  for  some  target  applications  it  enables  comparison  of  schemes  on 
an  even  footing. 

This  paper  is  organized  as  follows.  Section  2.1  provides  a  brief  overview  of  the 
Interacting  Multiple  Model  Kalman  Filter  used  here.  The  general  theory  of  the  IMMKF 
is  presented  in  [Bar].  A  more  detailed  description  of  this  filter  is  contained  in  [Kas96a]. 
Section  2.2  shows  how  the  expected  discrimination  gain  can  be  approximately  computed 
for  such  filters.  Section  2.3  presents  details  on  the  optimization  scheme  used  to  compute 
the  globally  optimum  assignment  of  sensors  to  targets.  Section  3  describes  a  test  scenario 
and  presents  simulation  results.  Section  4  presents  conclusions  and  suggestions  for 
further  work. 

2.1  Interacting  Multiple  Model  Kalman  Filters 

This  section  provides  a  brief  overview  of  Interacting  Multiple  Model  Kalman  Filters 
(IMMKF)  and  defines  notation  required  in  discrimination  gain  calculations.  The  IMMKF 
blends  motion  models  matched  to  different  flight  regimes  to  achieve  improved 
performance  against  maneuvering  targets  and  minimize  latency  due  to  model  switching. 
As  part  of  the  design  process  models  are  selected  to  approximate  different  flight  regimes 
for  the  target.  The  IMMKF  employs  a  soft  switching  scheme  that  probabilistically  blends 
the  output  of  the  different  filters.  The  probability  /n : .  ( k )  that  model  j  is  active  at  time 

k  is  estimated  as  part  of  the  IMMKF.  For  this  feasibility  study  we  use  an  IMMKF  based 
on  one  model  for  uniform  straight  and  level  flight  with  small  random  accelerations  and 
another  model  tuned  for  a  constant  rate  turn,  together  with  random  changes  in  the  turn 
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rate.  For  more  realistic  applications,  it  may  be  appropriate  to  include  additional  models 
for  extreme  acceleration  maneuvers.  Model  blending  is  achieved  assuming  that  the  target 
undergoes  transitions  between  flight  modes  modeled  by  a  Markov  chain.  The  estimates 
for  the  different  models  are  combined  with  the  residual  errors  for  the  models  to  estimate 
the  posterior  probability  that  each  mode  is  active. 

The  IMMKF  is  based  on  a  stochastic  hybrid  system  model  with  a  state  evolution  equation 
of  the  form 


X(k)  =  f[k  - 1,  X(k  - 1),  m(k)]  +  g[m(k),  k- 1,  X(k  - 1),  \(k,  m(k))\ , 


where  k  is  the  time  index;  X(k)  is  a  vector  describing  the  state  of  the  system  at  time  k 
(e.g.  position,  velocity,  turn-rate);  /  is  the  state  transition  function;  g  is  the  process  noise 
function;  m(k )  is  the  system  mode  during  the  interval  prior  to  time  k  (e.g.  uniform 
motion,  or  constant-rate  turn  modes);  \(k,  m(k))  is  a  process  noise  vector.  The  process 
noise  g  generally  depends  on  the  mode  m.  For  example,  in  straight  and  level  flight,  the 
process  noise  models  the  effect  of  small  inaccuracies  in  guidance  while  for  a  constant  rate 
turn,  the  process  noise  also  models  changes  in  the  turn  rate.  The  system  mode  m  is 
modeled  as  a  random  event  with  Markov  dynamics.  The  IMMKF  used  here  has  two 
modes:  mu  for  nearly  uniform  motion  and  m,  for  a  nearly  constant  rate  turn.  For 
example,  mu  (k  - 1)  denotes  the  event  where  a  target  is  undergoing  uniform  motion  at 
time  k  —  1  and  mT(k)  denotes  the  event  where  it  is  turning  at  time  k  . 

A  summary  of  the  IMMKF  and  its  notation  is  shown  in  Table  1  and  its  data  flow  is  shown 
in  Figure  2.  The  basic  operation  of  the  two-model  filter  comprised  of  a  uniform  motion 
filter  and  a  constant  turn  rate  filter  is  summarized  as  follows. 
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Table  1 

Interacting  Multiple  Model  Kalman  Filter  Summary 


Estimate  Mixing  (for  model  j ): 

model  transition  probability  for  time  interval  T :  Knj{T) 
probability  of  model  j  at  time  k  —  1 :  jUj(k-l) 
predicted  model  probability:  /r  =  ^  nj;H,  ( k  - 1) 

i 

model  mixing  probability:  junj  =  Tt^/u^k  -  1)  /  //  , 
mixed  state  estimate  for  model  j  : 

x0j(k  -  \\k  - 1)  =  £  x,(k -\\k- 1 )//,, 

i 

mixed  state  variance  estimate  for  model  j  : 

i 

Vtj  =  [x,.(k  -  II  k  - 1)  -  X0j(k  -II  k-  l)][x,.(k  -  \\k  - 1)  -  X0j(k  -II  k- 1)] T 


Filtering  (for  each  model  j ): 

state  prediction:  Xj(k\  k  - 1)  =  F j(k  -  l)X0y.(k  -II  k- 1) 
covariance  prediction: 

Pj(k\  k-  1)  =  F  j(k  -  1)P0  j(k\  k  -  1)F  .(k  -  1)T  +  G  .(k  -  1)Q  .(k  -  1)G  .(k  -  1)T 


measurement  residual:  Vj(k)  =  z (k)  -  H  j(k  -  1)X  j(k\  k  - 1) 

residual  covariance  estimate:  S j(k)  =  H  j(k  -  l)Py(kl  k  -  1)H  j(k  -  1)T  +  R ;.( k  - 1) 

filter  gain:  W )(k)  =  Pj(k\k  -  1)H  .(k  -  1)TS j(k)~l 

measurement  update:  Xj(k\k)  =  Xj(k\k  - 1)  +  Wj(k)Vj(k) 

covariance  update:  P.(fclfc)  =  Pj(k\k  -1)  -  W;(^S  ;(/:)W;(^)T 

Model  probability  update: 


model  likelihood  function:  A  =  (2 n) 


-1 


-1/2 

S/O')  exp 


-l-vAkfSAjr'vAk) 


model  probability:  ju j  =  // .  A .  /  ^  A,  A(. 


Estimate  combination: 

state  estimate:  X(kl  k)  =  ^  X;.(kl  /c )//; 

j 

covariance  estimate: 


P(fcl£)  =  I/^.  P^.(klk)  + 


X(k\k)-X.(k\k) 


T 


X(klk)-X^.(kl&) 


5 


1.  Estimate  Mixing:  To  initiate  IMMKF  processing  of  a  measurement,  the  previous  state 
estimates  for  the  models  are  blended  to  produce  mixed  state  estimates  (see  Table  1).  This 
results  in  a  modification  of  the  target  prediction  equations  and  accounts  for  mode 
switching  between  the  time  of  the  last  track  update  and  the  generation  of  the  current 
report.  For  example,  in  this  process  some  of  the  turning  model’s  estimate  is  included  in 
the  input  to  the  uniform  motion  filter,  with  appropriate  weighting.  This  accounts  for  the 
possibility  that  a  target  that  was  turning  at  the  conclusion  of  the  previous  measurement 
cycle  transitions  to  uniform  motion  for  the  current  measurement  cycle.  The  key  element 
in  this  computation  is  the  model  mixing  probability  junj .  This  is  the  probability  that  the 

target  was  in  mode  i  for  the  previous  scan,  given  that  it  is  now  in  mode  j  .  This  is  used 
to  construct  the  mixed  state  estimate  X0/  and  mixed  covariance  estimate  P0 .  for  each 
model  j  (see  Table  1).  In  the  mixed  covariance  estimate,  the  factor  IT  models  increased 
uncertainty  due  to  disagreement  between  the  model  estimates. 

2.  Filtering:  Each  mixed  state  estimate  is  updated  with  the  current  radar  report.  The 
standard  Kalman  filter  equations  are  used  with  the  appropriate  target  motion  models  and 
plant  noise  models.  This  produces  target  state  estimates  conditioned  on  the  target  model. 

3.  Model  Probability  Updating:  After  each  model  has  been  updated  with  the  radar  plot, 
new  model  probabilities  are  computed  from  the  model  likelihood’ s  A.  given  in 

Table  1  and  the  predicted  model  probability  //X  computed  during  the  preceding  mixing 
step. 

4.  Estimate  Combination:  The  state  and  covariance  estimates  for  all  of  the  models  are 
combined  using  the  updated  model  probability  weights  // . .  There  is  an  extra 

contribution  to  the  covariance  P(klk)  generated  by  disagreement  between  the  models 

[A  A  ir  A  A  "j  T 

X(k\k)  -  Xj(k\k)  X(klfc)  -  X;(klk)  ,  similar  in  form  to  the  mixed  state 

covariance  estimate  P0j(k\k).  Each  separate  model  state  is  maintained  internally  by  the 

system.  The  combined  estimate  is  only  required  for  output  of  the  system  track  to  other 
elements  of  the  air  defense  system. 

The  important  quantity  required  for  discrimination  gain  calculations  is  the  IMMKF 
combined  state  density  at  the  conclusion  of  each  measurement  cycle  given  by 
p(X,  j\Z)  =  jUjp(X\Z)  where 


p(XlZ)  =  \27iP{k\ k)\~in  exp 

V 


1 

2 


(X  -  X(k\ k))TP(k\k)~l(x  -  X(k\  k)) 


x 


with  the  other  quantities  defined  in  Table  1. 

2.2  Discrimination  Gain  and  the  IMMKF 
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To  compute  the  expected  discrimination  gain  using  an  IMMKF  filter,  suppose  that  we 
have  a  set  of  measurements  Z  made  up  to  time  tk_{ .  The  basic  idea  is  to  compute  the 
expected  discrimination  between  the  predicted  density  when  no  observation  is  made  and 
the  density  obtained  when  a  particular  sensor  is  used.  For  multivariate  Gaussian  qQ(X) 
and  q,(X)  with  means  X0  and  X,  and  variances  E0  and  Z, ,  the  discrimination  q1  with 
respect  to  q0  is  [Blahut] 


£(«,;<7ci)  =  |  <®9i(X)iog(^,(X)/^0(X)) 


e»+(x, 


Xq)(X,  -  x„)T) 


To  evaluate  the  relative  utility  of  an  observation  made  at  time  tk ,  first  note  that  if  no 
observation  is  made  then  the  target  state  density  is  obtained  from  the  IMMKF  prediction 
equations  alone.  Let  pQ(X,  j\ Z)  =  juJ p0(X\ Z)  denote  the  predicted  density.  Then 


f 

p0(X\ Z)  =  \2riP(k\  k  - 1)|"1/2  exp 

v 


1 

2 


(X  -  X(k\  k  -  l))TP(k\  k  - 1)"1  (X  -  X(k\  k  - 1)) 


x 

) 


where  P(k\k-l)  =  ^JjujPj(k\k-l) ,  X(k\k  - 1)  =  ^JUj  Xj(k\k  - 1)  and  is  given  in 

j  j 


Table  1. 


Now  if  a  new  measurement  z  obtained  with  a  particular  sensor  is  received  at  time  tk ,  the 
new  observation  set  is  Z'=j?}uZ.  The  IMMKF  computes  the  joint  conditional  density 
for  a  target  to  be  in  state  X  and  flight  mode  j ,  p(X,  j\Z )  =  jUjp(X\Z) .  The 
discrimination  of  p(X,  j\Z)  with  respect  to  p0(X,  j\Z)  is 

L  =  £  {  dXp(j,  XI Z')  In  (p(j,  XI Z')  /  Po(j,  XI Z)) 

j 

Note  that  this  depends  on  the  particular  observation  z  ■  For  sensor  management,  we 
require  the  expected  value  of  this  quantity  with  respect  to  the  conditional  density  p{z\Z) . 
This  is 

SL  =  Y,\ cIXclzp(z\Z)p(j,X\Z')\n{p(j,X\Z')/ Po(j,X\Z)) 

j 

The  integral  dz  over  the  new  observation  density  is  complicated.  However  p(z\Z)  is 
sharply  peaked  about  the  expected  value  of  the  observation  for  scan  k 
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zm  1  =  H^jUj  Xkjk  1  .  To  lowest  order  we  may  approximate  this  distribution  by  a  delta 

j 

function  at  zm~l .  Then  defining  Z'  =  jf“''juZ,  we  can  approximate 

£  L  »  £  jdX p(j,XIZ')ln(p(y, XIZ')  /  pQ(j,  XI Z)) . 

j 

Since  p(XIZ)  and  fi]  are  separately  normalized 

::L  =  ln(//;  / //0;)  +{ dXp(XIZ')ln(p(XIZ')/  P0(XIZ)) 


Using  the  fact  that  p(X,  j'IZ)  and  p0(X,  jlZ)  are  Gaussian,  we  obtain 


L  =  £{//,  ln(//;  / /i0;) 

j 

P(k\k-  lr^P(iklik)  -  P(klk  - 1)  +  (x(fclfc)  -  X(iklik  -l))(x(Ufc)  -  X(ifclik  - 1))' 


1 

H — tr 
2 


1,  P(fclfc) 

- log - 

2  P(k\k-\) 


The  gain  :  L  depends  on  which  sensor  is  used  through  the  sensor  covariance  dependence 
of  P(TI  k) .  This  expression,  evaluated  for  each  target  sensor  pair,  is  the  discrimination 
gain  used  here  as  a  sensor  effectiveness  metric. 


2.3  Optimization 

Given  the  gains  £  L  computed  above,  we  must  now  assign  sensors  to  targets.  The 
sensors  are  indexed  5  =  1, . . . ,  S  and  the  targets  are  indexed  t  =  1, . . . ,  T .  Each  sensor  has  a 
fixed  capacity  Ts .  This  is  the  maximum  number  of  targets  that  can  be  sensed  on  each 
sensor  scan.  The  discrimination  gain  when  sensor  s  is  assigned  to  target  t  is  denoted 
Gst.  Our  objective  is  to  maximize  the  total  gain  across  all  of  the  targets. 

The  optimization  problem  considered  in  this  paper  is  one  of  assigning  sensors  to  targets 
in  such  a  way  that  the  one  step  ahead  discrimination  gain  of  the  track  set  is  maximized. 

To  do  this,  the  one  step  ahead  predicted  discrimination  gain  G(-  -of  a  track  j  after  it  is 

updated  with  covariance  data  from  sensor  i  is  determined  for  all  pairs  i,j.  This  can  be 
formulated  as  a  Linear  Programming  problem  with  suitable  constraints.  The  next  step  in 
the  formulation  is  to  consider  the  situation  where  more  than  one  sensor  may  be  assigned 
to  the  same  target.  As  has  been  previously  pointed  out  [Nash],  one  method  is  to  construct 
pseudo  sensors  comprised  of  combinations  of  the  basic  sensors.  This  allows  any 


combination  of  sensors  to  be  assigned  to  a  single  target  simply  by  considering  that  a 
single  pseudo  sensor  has  been  assigned  to  it.  The  number  of  "sensors"  is  thus  increased 
from  S  to  2s  -  1.  The  problem  is  now  to  make  an  assignment  of  these  pseudo  sensors  to 
the  targets  in  an  optimal  way. 

There  are  however  constraints  on  this  assignment.  One  of  the  most  important  is  the 
maximum  tracking  capacity  of  a  sensor.  That  is,  given  a  specified  time  interval,  some 
sensors  can  scan  a  certain  volume  of  space  and  also  track  a  specified  number  of  targets. 
They  cannot  exceed  this  maximum  tracking  capacity  in  the  specified  time  period.  If  this 
maximum  tracking  capacity  is  known  for  each  of  the  basic  sensors,  then  this  must  be 
accounted  for  when  the  pseudo  sensors  are  assigned.  Previous  workers  have  not  handled 
this  accounting  correctly  [Nash]  when  they  define  additional  maximum  tracking 
capacities  for  the  pseudo  sensors.  Surprisingly, The  original  maximum  tracking  capacities 
are  all  that  is  needed. 

The  constraints  can  be  handled  exactly  as  follows.  Let  the  basic  sensors  be  numbered 

from  1  to  S.  Let  the  pseudo  sensors  be  numbered  from  S  +  1  up  to  2s  -  1.  For  each  basic 
sensor  k,  let  J(k)  be  the  set  of  integers  consisting  of  k  and  the  integer  numbers  of  the 
pseudo  sensors  which  contain  sensor  k  in  their  combination.  There  will  be  2  '  integers  in 
each  set  J(k).  These  sets  J(k)  will  appear  in  the  constraints  that  are  given  for  the  Linear 
Programming  formulation  to  our  optimal  assignment  problem. 

Example  1. 

For  example,  with  three  basic  sensors,  s  =  3,  and  2s  -1=7.  Let  S1,S2,S3,S4,S5,S6,  and 
S7  be  the  designations  of  the  sensors  with  S4  =  [S1,S2],  S5  =  [SI, S3},  S6  =  [S2,S3], 
and  S7  =  {S1,S2,S3}.  The  integer  sets  J(k)  then  contain  2s  4  integers  and  are  then:  J(l) 
=  { 1,4, 5, 7],  J(2)  =  [2, 4, 6, 7],  J(3)  =  [3, 5, 6, 7}. 


We  can  now  state  the  linear  programming  problem. 

2s  -1  t 

maximize  C  =  EE®.*, 


<=1  7=1 


subject  to  the  constraints 


(2.3.1) 


1  forj  =  l,  L  ,t  (2.3.2) 

i= 1 

t 

V  'Y'xij<Tk  fork  =  l,  L  ,s  (2.3.3) 

i  ej(k)  j= 1 

A'- j  >  0  for  all  pairs  ij  (2.3.4) 

and  where  Tk  is  the  maximum  tracking  capacity  of  the  basic  sensor  k.  In  the  LP  solution, 
each  xst  will  be  either  0  or  1 .  When  xst  =  1 ,  sensor  5  is  assigned  to  target  t . 
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Equation  (2.3.2)  fixes  each  j  and  sums  over  the  sensors  including  the  pseudo  sensors. 

The  fact  that  this  is  constrained  to  be  less  than  or  equal  to  one  means  that  only  one  or 
none  of  the  sensors  (recall,  the  pseudo  sensors  are  combinations)  will  be  assigned  to 
target  j.  Equation  (2.3.1)  clearly  sums  the  Gtj  for  which  there  is  a  sensor  assigned 

because  the  x(  ;  are  zero  or  one.  Equation  (2.3.3),  by  virtue  of  the  sets  J(k),  insures  that 

the  number  of  targets  tracked  by  the  assignment  sensors  does  not  exceed  the  maximum 
tracking  capacity  of  the  basic  sensor  k.  If  there  are  fewer  targets  than  the  maximum 
sensor  tracking  capacity,  then  some  pseudo  sensors  will  be  assigned  to  individual  targets. 
If  there  are  more  targets  than  the  total  tracking  capacity,  then  Equation  (2.3.2)  allows 
some  targets  to  not  have  sensors  assigned. 

3.  Test  Results  Typical  test  results  for  a  2  sensor,  8  target  test  case  are  shown  in  Figure 
1-3.  For  this  test  case,  the  two  sensors  are  at  fixed  positions.  Sensor  locations  and  other 
system  parameters  are  shown  in  Table  2.  Two  scenarios  were  examined.  The  first  test 
simulates  the  situation  where  two  sensors  must  track  8  closing  targets.  For  the  second 
case,  the  sensor  geometry  is  changed  to  reflect  targets  transiting  along  a  political 
boundary.  The  sensors  provides  measurements  to  the  central  fusion  center  at  1  second 
intervals,  referred  to  as  a  scan.  In  each  scan,  we  assume  that  the  sensor  can  observe  two 
targets,  so  half  the  targets  are  observed  on  each  scan.  The  targets  are  numbered  1-8,  with 
target  1  initially  the  southern-most  target  and  target  8  the  northern-most.  Target  4 
undergoes  a  standard  rate  turn.  The  remaining  targets  are  in  straight  line  flight.  Figure  3 
shows  the  true  turn  rate  and  the  turn  rate  estimate  for  Target  4  averaged  over  10  runs  of 
the  scenario.  There  is  some  lag  in  the  filter  and  it  correctly  estimates  the  peak  turn  rate  at 
about  .03  radians/sec. 

Figure  3  shows  the  expected  discrimination  gain  vs.  scan  k  averaged  over  the  same  10 
runs  of  the  test  scenario.  The  upper  curve  is  the  expected  gain  of  the  optimal  sensor  for 
the  maneuvering  target.  The  lower  curve  is  the  optimal  gain  averaged  over  the  7  non¬ 
maneuvering  targets.  As  expected,  the  maneuvering  target  generally  has  greater 
uncertainty  and  has  greater  expected  discrimination  gain  as  a  result.  This  is  confirmed  by 
the  fact  that  the  maneuvering  target  is  observed  nearly  60%  of  the  time,  while  the  other 
targets  are  observed  somewhat  less  than  50%  of  the  time.  Note  that  even  though  there  is 
significant  turn-rate  lag  as  shown  in  Figure  2,  the  expected  gain  begins  to  grow  as  soon  as 
the  target  initiates  its  maneuver  in  scan  20. 

To  quantify  the  efficacy  of  discrimination  gain  as  a  measure  of  sensor  effectiveness  we 
computed  the  RMS  position  error  and  total  error  covariance  for  three  alternative  sensor 
management  schemes.  These  are  given  in  Tables  3  and  4.  The  rows  labeled  "DG"  use 
discrimination  gain  with  the  simplex  optimization.  The  rows  labeled  "Fixed"  use  a  fixed 
round-robin  sensing  schedule.  There  are  many  possible  alternatives.  We  chose  to  use  the 
schedule  where  sensor  1  observes  target  pair  (1,5),  then  (2,6),  (3,7  )  while  sensor  2 
observes  targets  (3,7),  then  (4,8)  and  so  on.  This  is  simple  to  implement  and  spreads  the 
observations  out  among  the  targets.  This  coordinates  sensing  but  does  not  adaptively 
compute  the  sensor  effectiveness.  For  the  rows  labeled  "Rand"  each  sensor  selects  targets 


10 


at  random  for  each  scan  cycle.  This  is  the  worst  case  approach  in  that  it  assumes  that  the 
sensing  is  completely  uncoordinated. 

In  order  to  establish  a  baseline  for  comparison,  we  assumed  that  all  of  the  scheduling 
schemes  maintain  central  tracks  updated  via  report  fusion  as  soon  as  measurements  are 
available.  Although  this  may  be  unrealistic  for  some  applications,  it  enables  simple  direct 
comparison  of  the  three  schemes. 

We  examined  two  measures  of  effectiveness  for  this  sensor  management  study.  The  first 
is  simply  the  RMS  position  error  for  each  target  after  each  measurement  update.  For 
scans  where  a  particular  target  is  not  updated,  its  predicted  position  X(k\  k  - 1)  is 
computed  and  compared  with  ground  truth  to  obtain  the  RMS  error. 

For  some  tracking  applications  other  quantities  such  as  the  velocity  accuracy  are  of 
primary  interest.  A  simple  metric  that  is  also  sensitive  to  this  more  global  aspect  of  the 
problem  can  be  obtained  from  the  covariance 

C  =  (X(k)  -  X(k\  k)){x(k)  -  X(k\  k)f , 

which  is  a  5x5  matrix.  The  determinant  has  units  [|C|]  =  in  /  s'  and  is  proportional  to 
the  volume  of  the  uncertainty  ellipsoid  of  the  target  estimate.  The  volume  of  total 
uncertainty  ellipsoid  for  the  8-target  problem  (in  a  40-dimensional  space)  is  proportional 
to  the  product  of  the  single  target  volumes. 

Global  discrimination  optimization  yields  the  best  performance  with  80  m  RMS  position 
error  averaged  across  all  the  targets  while  random  assignment  yields  the  worst  with  90  m 
average  RMS  error.  Also,  the  greatest  gain  was  obtained  against  the  maneuvering  target 
where  discrimination  yielded  93  m  RMS  position  error  while  the  random  approach 
suffered  from  a  109  m  RMS  error  for  the  maneuvering  target.  Similar  results  were 
obtained  for  both  test  scenarios.  These  qualitative  results  are  reflected  in  the  covariance 
determinants  for  these  scheduling  methods,  as  well.  Interestingly,  the  round-robin  sensor 
scheduling  was  able  to  provide  significant  performance  improvement  relative  to  the 
random  scheduling,  although  its  performance  still  fell  below  that  of  discrimination  gain. 
This  serves  to  emphasize  the  utility  of  multisensor  scheduling. 
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Figure  1  -  XY  Plot  (meters)  of  test  scenario  for  a  collection  of  targets  closing  with  a  pair 
of  sensor  platforms  Si  and  S2  showing  target  trajectories  and  sensor  detections.  All 
targets  enter  from  the  left.  Target  4  (indicated)  executes  a  standard  rate  turn  during  the 
scenario. 


CO 


(rad/sec) 


Figure  2  -  IMMKF  Turn  Rate  Estimate  Turn  rate  (dashed)  and  turn  rate  estimate  for 
target  4  averaged  over  10  scenario  runs. 
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Figure  3  -  Expected  Discrimination  Gain  with  Global  Optimization  Expected 
discrimination  is  greater  for  maneuvering  target  leading  to  more  sensor  resource 
allocation  and  improved  tracking  performance.  Expected  discrimination  for  maneuvering 
target  (target  4)  -  upper  curve  (solid  line).  Average  expected  discrimination  for  non¬ 
maneuvering  targets  -  lower  curve  (dashed  line).  Both  curves  are  average  results 
obtained  over  10  runs  of  the  closing  target  scenario.  Note  that  the  gain  for  the 
maneuvering  target  increases  relative  to  the  non-maneuvering  targets  as  soon  as  the 
maneuver  is  initiated  at  scan  20. 
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ag  =  1°  Sensor  1  location  (x,y) 

1 1  km  -5  km 

Sensor  2  location  -  closing  targets  (x,y) 

11  km  11  km 

Sensor  2  location  -  transiting  targets  (x,y) 

-5  km  -5  km 

Sensor  update  rate 

1  sec'1 

Target  speed 

100  m/sec 

<7 

range 

100  m 

1° 

Table  2  -  Sensor  and  target  parameters  used  to  test  sensor  effectiveness  measure. 


Tgtl 

Tgt2 

Tgt3 

Tgt4 

Tgt5 

Tgt6 

Tgt7 

Tgt8 

Avg 

DG  RMS  (m) 

85 

76 

83 

93 

76 

74 

73 

75 

80 

DG  Det 

25 

14 

19 

1470 

14 

14 

454 

10 

252 

Fixed  RMS 

90 

76 

81 

110 

78 

76 

79 

84 

85 

Fixed  Det 

32 

14 

19 

2824 

19 

14 

525 

25 

434 

Rand  RMS 

82 

90 

83 

109 

83 

95 

95 

80 

90 

Rand  Det 

19 

19 

19 

2293 

19 

32 

902 

19 

415 

Table  3  —  Test  results  obtained  for  10  runs  of  2-sensor  closing  target  scenario.  Use  of 
discrimination  gain  (DG)  to  select  the  optimal  sensor  assignment  reduces  the  average 
track  error  as  measured  by  both  the  RMS  position  error  and  the  determinant  of  the 
average  covariance  (highlighted  column  at  right).  DG  is  compared  to  a  fixed  sensor 
schedule  and  to  a  random  sensor  schedule  (Rand)  that  simulates  the  situation  where  there 
is  no  coordination  between  the  sensors.  DG  provides  the  most  benefit  against  target  4 
which  maneuvers  during  the  scenario. 


Tgtl 

Tgt2 

Tgt3 

Tgt4 

Tgt5 

Tgt6 

Tgt7 

Tgt8 

Avg 

DG  RMS  (m) 

77 

81 

80 

102 

87 

73 

80 

88 

84 

DG  Det 

11 

20 

14 

1850 

25 

7.6 

454 

24 

300 

Fixed  RMS 

91 

86 

85 

111 

78 

77 

85 

91 

89 

Fixed  Det 

32 

19 

19 

2824 

14 

11 

19 

32 

371 

Rand  RMS 

85 

94 

92 

116 

86 

105 

93 

90 

95 

Rand  Det 

19 

32 

25 

2824 

19 

52 

792 

19 

472 

Table  4  -  Test  results  obtained  for  10  runs  of  transiting  target  scenario  are  qualitatively 
similar  to  results  for  closing  target  scenario  (Table  3). 
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4.  Discussion  To  summarize,  we  have  demonstrated  that  the  formalism  of  the  IMMKF 
can  be  used  as  a  framework  for  computing  the  expected  information  gain  when  multiple 
sensors  observe  a  collection  of  airborne  targets.  When  combined  with  the  IMMKF, 
discrimination  gain  is  able  to  recognize  target  maneuvers  and  respond  by  allocating 
additional  sensor  resource  to  the  maneuvering  target.  The  sensor  target  assignment  that 
optimizes  the  total  discrimination  across  all  of  the  targets  in  a  surveillance  volume  can  be 
readily  computed.  Application  of  this  optimal  assignment  results  in  improved  tracking 
performance  relative  to  random  assignment  of  sensors  to  targets  or  fixed  round-robin 
scheduling. 

The  IMMKF  used  in  this  assessment  was  designed  for  an  air  traffic  control  application. 
For  use  in  air-to-air  tracking  applications,  different  versions  of  the  IMMKF  must  be 
developed  to  treat  target  motion  in  3-dimensions  and  to  incorporate  range-rate 
measurements.  Another  topic  that  must  be  examined  is  the  use  of  simpler  filters  that  can 
still  support  discrimination  gain  evaluation.  For  example,  the  adaptive  single  model 
Kalman  filter  can  also  be  used  in  this  application  but  was  not  studied  here.  It  may  be  that 
for  some  applications  such  as  tracking  low-priority  targets,  the  computational  simplicity 
of  the  single  model  filter  outweighs  its  disadvantages  relative  to  the  IMMKF.  This  study 
has  assumed  that  the  target  detection  probability  is  1  and  that  there  is  no  clutter. 
Furthermore,  it  was  assumed  that  the  track  is  a  confirmed  target.  However,  if  the  track  is 
a  tentative  target,  then  additional  information  is  obtained  when  the  target  is  observed  with 
the  sensor.  These  effects  can  be  included  in  the  discrimination  gain  calculations. 

One  might  expect  the  improved  performance  for  the  maneuvering  target  to  come  at  the 
expense  of  some  loss  of  performance  for  the  non-maneuvering  targets.  However,  this  is 
not  the  case.  By  improving  sensor  allocation  across  all  of  the  targets,  average 
performance  is  improved  for  both  the  maneuvering  and  non-maneuvering  targets. 

Finally,  this  study  has  addressed  sensor  effectiveness  independent  of  the  tactical  or 
strategic  utility  of  knowledge  about  the  targets.  In  order  to  address  this  aspect  of  the 
problem,  discrimination  must  be  combined  with  a  means  of  assessing  preferences  as  well. 
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